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Physical Nucleon Properties from Lattice QCD 
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We demonstrate that the extremely accurate lattice QCD data for the mass of the nucleon recently 
obtained by CP-PACS, combined with modern chiral extrapolation techniques, leads to a value for 
the mass of the physical nucleon which has a systematic error of less than one percent. 
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Hadronic physics presents fascinating theoretical chal- 
lenges to the understanding of strongly interacting sys- 
tems in terms of their fundamental degrees of freedom in 
QCD, quarks and gluons. Traditional perturbative meth- 
ods in quantum field theory cannot account for the highly 
non-trivial structure of the QCD vacuum - which is dom- 
inated by large scale quark and gluon correlations. Lat- 
tice gauge theory I'l has so far provided the only rigorous 
method of solving non-perturbative QCD. We will show 
that recent progress within the field , together with ad- 
vances in effectively field theory (EFT) now permit 
the accurate, model-independent extraction of the phys- 
ical properties of hadrons from lattice QCD simulations, 
even though it is not yet feasible to make calculations at 
the physical quark mass. 

It is now possible to make extremely accurate calcula- 
tions of the masses of hadrons, such as the nucleon, with 
dynamical fermions. Indeed the CP-PACS group has just 
reported data with a precision of order 1% However, 
such precise evaluations are limited to quark masses an 
order of magnitude larger than those found in nature. 
In order to compare with experiment, which is after all 
one of the main aims in the field, it is therefore neces- 
sary to extrapolate in quark mass. It is well known that 
such an extrapolation is complicated by the unavoidable 
non-analytic behaviour in quark mass, which arises from 
Goldstone boson loops in QCD with dynamically broken 
chiral symmetry 

Early work motivated by the important role of Gold- 
stone bosons led to the construction of chiral quark mod- 
els [3, which incorporated this non-analytic behavior. 
An alternative, systematic approach, designed to avoid 
reference to a model, involved the construction of an ef- 
fective field theory (EFT) to describe QCD at low en- 
ergy 0. The application to baryons has developed to 
the point where chiral perturbation theory (xPT) is now 
understood as a rigorous approach near the chiral limit 

Hi. 

Because it is defined as an expansion in momenta and 
masses about the chiral limit, xPT provides an attrac- 
tive approach to the problem of quark mass extrapola- 
tion for lattice QCD. The advantages of formulating xPT 
with a physical regularisation scheme, as opposed to the 



commonly implemented dimensional regularisation, have 
been demonstrated by Donoghue et al. 9]. Early imple- 
mentations of a finite-range regulator (ERR) to evaluate 
chiral loop integrals in xPT suggested that, in the con- 
text of the extrapolation of lattice data from relatively 
large quark masses, ERR provides a more reliable proce- 
dure [l(j| . There has been considerable debate on whether 
current lattice data is within the scope of dimensionally- 
regularised xPT or whether the form of the ERR chosen 
introduces significant model-dependence 0|. However, 
this issue has now been resolved through a recent de- 
tailed study of numerous regularisation schemes in xPT, 
both dimensional and ERR. This study quantified the 
applicable range of the EET 3] and established that all 
of the ERR considered provided equivalent results over 
the range ml < 0.8 GeV^. 

Here we demonstrate that by adopting the ERR formu- 
lation of xPT one can improve its convergence properties 
to the point where the chiral extrapolation to the physi- 
cal quark mass can be carried out in a model-independent 
way with a systematic uncertainty of less than one per- 
cent. Given this remarkable result, the current errors 
on extrapolated quantities are dominated by the sta- 
tistical errors arising from the large extrapolation dis- 
tance - the lightest simulated pion mass being typically 
ml ~ 0.27 GeV^, in comparison with the physical value, 
0.02 GeV^. 

In the usual formulation of effective field theory the 
nucleon mass as a function of the pion mass (given that 
ml oc TOg 

M 



^3|) has the formal expansion: 
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In principle the coefficients, a„, can be expressed in terms 
of the parameters of the underlying effective Lagrangian 
to a given order of chiral perturbation theory. In prac- 
tice, for current applications to lattice QCD, the param- 
eters must be determined by fitting to the lattice results 
themselves. The additional terms, ctntt and ctait, are 
the two self-energy terms involving the (Goldstone) pion 
which yield the leading and next-to-leading non-analytic 
behaviour of M^. As these terms involve the coupling 
constants in the chiral limit, which arc essentially model- 
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independent the only additional complication they 
add is that the ultra-violet behaviour of the loop integrals 
must be regulated in some way. 

Traditionally, one uses dimensional regularisation 
which (after infinite renormalisation of all even powers 
of m^) leaves only the non- analytic terms - clna^^ and 
CNLNAnT'^ ln(m7r//x), respectively. Within dimensional 
regularisation one then arrives at a truncated power se- 
ries for the chiral expansion, 

Mn = Co + C2ml + CLNA"!^ -I- C4m^ 

+cnlnato^ hi — + CQml + . .. , (2) 
A* 

where the bare parameters, a^, have been replaced by 
the finite, renormalised coefficients, q. Through the chi- 
ral logarithm one has an additional mass scale, fj., but 
the dependence on this is eliminated by matching C4 to 
"data" (in this case lattice QCD). Provided the series 
expansion in Eq. (0) is convergent over the range of val- 
ues of TTT-Tr where the lattice data exists, one can then use 
Eq. (O to evaluate Mjv at the physical pion mass. Un- 
fortunately there is considerable evidence that this series 
is not sufficiently convergent [sl llllll^ Hslllil . 

The physical reason for the lack of convergence seems 
to be that, apart from the usual symmetry breaking scale 
in xPT, A^sB ~ 47r/^ ~ 1 GeV, there is another mass 
scale in the problem which does not appear explicitly 
in a dimensionally regularised theory. This additional 
scale. A, is the inverse of the size of the nucleon, the 
source of the pion cloud 0. Of course, this scale must 
appear implicitly in the correct effective field theory as 
a correlation between the coefficients c„. However, since 
we want to access pion masses above 0.5 GeV and the 
nucleon size of 0.5-1.0 fm suggests A ~ 0.4 GeV, it is 
clear that ttItt/A > 1 and the naive series should not be 
expected to converge over the range of pion masses of 
interest. Thisproblem can be overcome through the use 
of aFRR [101. 

Clearly, in extracting information from lattice data for 
physical hadrons one must avoid model dependence. In 
line with the implicit /^-dependence of the coefficients in 
the familiar dimensionally regulated xPT, in the case of 
a FRR the systematic expansion of the nucleon mass is: 

Mm = a^ + a^ml + a^ml 

-|-CTN7r(m^, A) + (TAx(m^, A), (3) 

where the dependence on the shape of the regulator is 
implicit. As for the mass scale, fj,, entering the usual 
chiral logarithms, the dependence on the value of A and 
the choice of regulator is once again eliminated, to the 
order of the series expansion, by fitting the coefficients, 
a^, to lattice data. The clear indication of success in 
eliminating model dependence and hence having found a 
suitable regularisation method, is that the higher order 
coefficients (for the best fit value of A) should be small 



and that the result of the extrapolation should be insen- 
sitive to the shape of the ultraviolet regulator. 

In order to investigate the model dependence associ- 
ated with the self-energies, several regulators are consid- 
ered. We evaluate the loop integrals in the heavy baryon 
limit 

167r^/^ Jo w(fc)(wBB' + w(fc)) 

taking u{k) to be either a sharp cut-off, 6{h. — k), a 
dipole, (1 -t- P a monopole, (1 + P / or fi- 
nally a Gaussian, exp(— fc^/A^). We note that the dipole 
is the most physical, in that the observed axial form fac- 
tor of the nucleon has this shape In Eq. Q for 
B = B' = N we have Gnn = 9\ (with 5^ = 1.26), 
while Gata ~ 32(7^/25, the SU(6) value. In addition, 
ijj{k) = ^k"^ + ml, UJNN = and ujna = 292 McV, the 
physical A-A^ mass splitting. Clearly the four choices of 
regulator have very different shapes, with the only com- 
mon feature being that they suppress the integrand for 
momenta greater than A. 

In addition to the FRR expansions of Eq. Q and the 
standard dimensionally regulated expansion of Eq. 
we also consider the case where Eq. (jJl is modified to 
maintain the square-root branch-cut at TOtt = i^AfA 
One can now compare the expansion about the chiral 
limit for these six different regularisation schemes in or- 
der to assess their rate of convergence. It turns out 
that all the FRR expansions precisely describe the di- 
mensional regularization expansion over the range to^ S 
(0,0.7) GeV^. Furthermore, the smooth, FRR formula- 
tions are consistent with each other, for the renormal- 
ized chiral coefficients, Co,2,4, to an extraordinarily pre- 
cise level . This ensures a reliable extrapolation to the 
regime of physical quark masses. 

We now exploit the accuracy of the most recent CP- 
PACS data H by using it to determine either the first 
four unknown parameters, co — ce, appearing in Eq. ^ 
(i.e. using dimensional regularisation, with and with- 
out the A — > Nn branch point) or the four parameters 
a.0,0,2, 0,4 and A in Eq. for each choice of ultra-violet 
regulator It is only the extreme accuracy of the 

data which makes the determination of as many as four 
parameters possible. Figure shows the resulting fits to 
the lattice data over the range G (0, 1.0) GeV^, with 
the corresponding parameters given in Tabled It is re- 
markable that three of the four curves based on Eq. ((SJ, 
the dipole, monopole and Gaussian cases, are indistin- 
guishable on this plot. Furthermore, we see from TablelH 
that the coefficient of in all of those cases is quite 
small - an order of magnitude smaller than the higher 
order coefficients obtained using Eq. which also os- 
cillate in sign. This indicates that the residual series, in- 
volving ao,2,4, is converging when the chiral non-analytic 
behaviour is evaluated with a FRR. 
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FIG. 1: Fits to lattice data for various ultra-violet regulators. 
The monopole, dipole and Gaussian cases are depicted by 
solid lines, indistinguishable on this plot, the sharp cutoff 
is shown by the dashed curve. The dimensional regularised 
forms are illustrated by the dash-dot curves, with the correct 
branch point corresponding to the higher curve. Lattice data 
is from Ref. igj. 



Regulator 


ao 


ai 


a4 




A 


Dim. Reg. 


0.700 


5.10 


3.87 


-2.35 




Dim. Reg. (BP) 


0.691 


5.32 


7.65 


-1.36 




Sharp Cutoff 


1.06 


1.06 


-0.249 




0.440 


Monopole 


1.38 


0.950 


-0.217 




0.443 


Dipole 


1.15 


0.998 


-0.227 




0.732 


Gaussian 


1.11 


1.02 


-0.234 




0.593 



TABLE L Bare, unrenormalised, parameters extracted from 
the fits to lattice data displayed in Figure (All quantities 
are in units of appropriate powers of GeV and = 1 GeV in 
Eq. ©.) 



As explained by Donoghue et al. p| , one can combine 
the order to°'^''* terms from the self-energies with the 
"bare" expansion parameters, ao,2,4j to obtain physically 
meaningful renormalized coefficients. These are shown 
in Table UTI in comparison with the corresponding coeffi- 
cients found using Eq. Q . Details of this renormalization 
procedure are given in Ref. Q ■ The degree of consistency 
between the best-fit values found using all three smooth 
choices of FRR is remarkably good, whereas the large nu- 
merical values and the change of sign in the coefficients 
found from the fit using Eq. suggest that the latter 
is nowhere near convergence. We can understand the 
problem very simply; it is not possible to accurately re- 
produce the necessary 1 /m^ behaviour of the chiral loops 
(for TOtt > A) with a 3rd order polynomial in m^. The 
final point to notice about the finite range regulators is 
that all but the sharp cutoff naturally yield a correction 
to the Goldberger-Treiman relation of approximately the 



Regulator 


Co 


C2 


C4 


Dim. Reg. 


0.700 


5.10 


3.87 


Dim. Reg. (BP) 0.750 


4.59 


6.43 


Sharp cutoff 


0.892 


3.16 


12.9 


Monopole 


0.914 


2.63 


26.1 


Dipole 


0.910 


2.72 


23.4 


Gaussian 


0.906 


2.79 


21.4 



TABLE IL Renormalised expansion coefficients in the chiral 
limit obtained from various regulator fits to lattice data. (All 
quantities are in units of appropriate powers of GeV.) 



right size. As shown by Birse and McGovern ,20], this 
correction dominates the non-analytic term of order 
which only arises at two loops in dimensionally regulated 
XPT. The absence of the correct term in the case 
of the sharp cutoff explains its slight deviation from the 
results of the other FRR. 

It is clear that the use of an EFT with FRR enables 
one to make a model-independent extrapolation of the 
nucleon mass as a function of the quark mass. Although 
minimal deviation is seen between the best fit curves, 
we need to determine how well these curves are in fact 
constrained by the statistical uncertainties of the lattice 
data. As all data points are statistically independent, the 
one-sigma deviation from the best-fit curve is defined by 
the region for (x^ — Xmin)/'^^/ < 1. We use a standard 
measure, weighted by the squared error of the simulated 
data point and Xmin corresponds to the optimum fit to 
the data. To preserve the infrared physics of EFT, we 
put a lower bound of 300 MeV (just above the N — A 
mass splitting) on the regulator parameter, A. 

We show the one-sigma variation from the best-fit 
dipole curve by the shaded region in Figure [3 The ex- 
trapolated values for the nucleon mass are shown in Ta- 
ble ^O] It is clear that the large lower-deviation on the 
extrapolated values are a result of the regulator mass 
not being constrained from above. Alternatively, for the 
FRR cases, one could fix A and vary an additional pa- 
rameter, og, in the spirit of dimensional regularisation. 
In that case the error bounds are symmetric and similar 
in magnitude to those found for the dimensionally regu- 
lated cases reported in Table ITlll The primary source of 
the large error band is the large extrapolation distance. 

We also quote the cr-commutator in Table IIIII - cal- 
culated as rn^dMN/drri^ at the physical pion mass [2ll |. 
Again the statistical errors are large. In this well 
as for the nucleon mass itself, it is clear that a simulated 
point at a pion mass ~ 0.1 GeV^ would greatly re- 
duce the statistical error in the extrapolation. We note 
that our study has focused on the systematic errors asso- 
ciated with chiral extrapolation, while those arising from 
finite lattice spacing and volume [isl . [2^ remain to be 
quantified. 
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the issue of model independence resolved, there is an ur- 
gent need for high precision lattice QCD simulations at 
^0.1 GeV^ in order to reduce the present statistical 
error on the extrapolation. Such simulations should be 
feasible with the new generation of lattice QCD comput- 
ers currently under construction. 

We thank M. Birse, J. McGovern and S.V. Wright for 
helpful conversations. This work was supported by the 
Australian Research Council. 



FIG. 2: Error analysis for the extraction of the nucleon mass 
using a dipole regulator. The shaded region corresponds to 
the region allowed within the present statistical errors. 



Regulator 


mN 




(TJV 


Dim. Reg. 


0.782l»;;g 


70+15 


Dim. Reg. (BP) 


U.OZO_o.i22 


65 


+ 15 
-15 


Sharp cutoff 


0.939l°:°I« 


41 


+35 
-12 


Monopole 


0.953l°:°f, 


35 


+36 
-8 


Dipole 


0.95lt»:°?^ 


36 


+35 
-13 


Gaussian 


0.9481°:°?° 


37 


+35 
-13 



TABLE III: Here we show the nucleon mass, mjv (GeV), and 
the sigma commutator, ajv (MeV), extrapolated to the phys- 
ical pion mass. 



To summarize, we have shown that the extremely pre- 
cise dynamical simulation data from CP-PACS permits 
us to determine four parameters in the chiral extrapo- 
lation formulas, Eqs. l(2Jl and Q. Whereas the former 
(involving dimensional-regularization) does not appear 
to be convergent over the required mass range, the im- 
proved convergence properties of a finite range regulator 
yield an excellent description of the data over the full 
mass range, regardless of the functional form chosen for 
the vertex regulator. Table Hill summarises the resulting 
values of the nucleon mass and the cr-commutator eval- 
uated at the physical pion mass. The nucleon mass is 
constrained by the chiral extrapolation technique based 
upon a FRR procedure to a systematic error less than 
1%. The cr-commutator is similarly constrained to within 
a systematic error of 3%. The values obtained are both 
consistent with earlier extractions from lattice data using 
a similar technique |2lj | and, within the relatively large 
statistical errors, in agreement with the experimental val- 
ues [13 . We note also that the systematic uncertainty in 
the extraction of the low energy constant, cq, is less than 
1%, while for C2 it is at the level of a few percent. With 
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